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ABSTRACT 

We obtained high angular resolution (~2") images of the 13 C0(J=1— >0) line 
and 2.7 millimeter continuum emission, and slightly lower resolution images of 
12 C0(J=1— >0) and C 18 0(J=1— >0) line emission toward the Herbig Ae star AB 
Aurigae. We resolve a circumstellar disk of diameter 780 AU (FWHM) with a 
velocity pattern consistent with a purely rotational disk at inclination 21.5° and 
position angle 58.6°. Using Keplerian disk models, we find a central source dy- 
namical mass of 2.8±0.1 M and a cutoff radius of 615 AU for the 13 CO emission. 
Inclination, mass, and radius determined from 12 CO and C 18 observations agree 
with those values, given optical depth and abundance effects. As a result of the 
high angular resolution of our observations, we confirm the existence of spiral 
structure suggested by near-IR scattered light images and show that the spiral 
arms represent density contrasts in the disk. 

Subject headings: circumstellar matter — stars: individual (AB Aur) 



1. Introduction 

Circumstellar disks often surround low (<2 M Q ) and intermediate (~2-10 M ) mass pre- 
main sequence objects, the T Tauri (TT) and Herbig Ae/Be (HAeBe) stars. Disk sizes are 
of order a few hundred AU with masses a few tenths of a solar mass or less (e.g. Beckwith & 
Sargent 1996; Natta, Grinin, & Mannings 2000). These values are similar to the proto-solar 
system (Weidenschilling 1977), suggesting that they may be the sites of planet formation. 
Disk temperature and density profiles, key properties for understanding how planets might 
emerge, have been inferred from spatially unresolved observations that rely on spectral en- 
ergy distributions (SEDs) and require assumptions about disk morphology (e.g. Kenyon & 
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Hartmann 1987; Dullemond, Dominik, & Natta 2001). However, higher angular and spectral 
resolution measurements of the gas and dust in disks are critical to quantifying the profiles 
directly and providing a context in which to interpret unresolved observations. Spatially 
and kinematically resolved images enable measurement of stellar mass, disk mass, radius, 
inclination (i), position angle (PA), and substructure (e.g. Koerner, Sargent, & Beckwith 
1993; Dutrey et al. 1998; Simon, Dutrey & Guilloteau 2001). Such observations of multiple 
spectral lines in DM Tau have allowed the exploration of vertical and radial disk structure 
(Dartois, Dutrey, & Guilloteau 2003). 

The HAe star AB Aurigae at 144 pc (Perryman et al. 1997) has been studied by nu- 
merous authors at optical, IR and radio wavelengths. Most recently, Semenov et al. (2004, 
hereafter S04) based a detailed analysis on millimeter observations at 5-12" resolution. While 
this angular resolution is insufficient to determine a stellar mass, their model fits lead to a 
disk radius of 400±200 AU, % of 17l|°, and PA of 80±30°. This inclination agrees well 
with estimates of i<20-30° from near-IR interferometry and scattered light (Millan-Gabet, 
Schloerb, & Traub 2001; Eisner et al. 2003, 2004; Grady et al. 1999; Fukagawa et al. 2004, 
hereafter F04). Earlier millimeter interferometry and recent mid-IR observations at angular 
resolution ~0'/5 suggest i~45-70° (Marsh et al. 1995; Mannings & Sargent 1997, hereafter 
MS97; Liu et al. 2004). 

High resolution interferometric images of circumstellar gas can determine the inclination 
definitively. It is likely that early studies of the AB Aur disk at moderate angular resolution 
were affected by spiral structure in the disk surface (F04) or large scale (~1000 AU) envelope 
emission (e.g. S04). Here we present high resolution (2'.'2) 2.7 mm continuum and 13 CO(1-0) 
line images and slightly lower resolution maps of the 12 CO(1-0) and C 18 O(l-0) emission from 
AB Aur. In addition to constraining the disk properties and central mass, these enable us 
to detect and probe the spiral features seen in the images of F04. 

2. Observations and Results 

Observations of AB Aur were carried out between January and May of 2004 at the 
Owens Valley Radio Observatory (OVRO) millimeter-wave array. Six 10.4 m antennas with 
cryogenically cooled SIS receivers were used in five array configurations with baselines ranging 
from 15 to 480 m. We configured the correlator to observe the emission lines of 13 CO at 
110.2 GHz and 12 CO at 115.3 GHz in 64 Hanning-smoothed channels of width 125 kHz, and 
the C 18 emission line at 109.8 GHz in 32, 250 kHz wide channels. The 12 CO and C 18 
observations have maximum baselines of 103 and 120 m respectively. Simultaneously, we 
obtained 2.7 mm continuum measurements over 7 GHz of bandwidth with OVRO's COBRA 
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correlator (Hawkins et al. 2004). The quasars J0530+135 and 3C111 were observed at 20 
minute intervals for phase and amplitude calibration. Absolute flux calibration, accurate 
to 10%, was based on measurements of Uranus and Neptune. Data reduction was carried 
out with MIR, an IDL-based package developed for the OVRO array by N. Scoville and J. 
Carpenter. Maps were made using MIRIAD (Sault, Teuben, & Wright 1995). 

For AB Aur, Figure 1 shows contours of integrated intensity for 12 CO (left panel), 
13 CO (center panel), and C 18 (right panel) emission at resolutions 4'/25, 3'/25, and 3'.'9, 
respectively overlaid on intensity- weighted mean velocity maps (color). Integrated line fluxes 
for 12 CO, 13 CO, and C 18 are 22.1, 3.3 and 0.78 Jy km s _1 over the velocity ranges 7.67 to 
4.27, 7.51 to 4.25, and 7.37 to 4.63 km s _1 , respectively. Outside these ranges, no emission 
was detected above the 3a level (360, 105, and 72 mJy/beam for 12 CO, 13 CO, and C 18 
respectively). We combined our 13 CO dataset with that of MS97 (re- reduced using MIR) 
to provide greater sensitivity and found the 13 CO flux agrees with that of MS97, within the 
uncertainties. To emphasize global spatial and velocity structure, the images were restored 
with circular beams equal in area to the naturally weighted beams. It is immediately evident 
that the 12 CO and 13 CO emission is resolved and the velocity varies smoothly across the maps. 
Due to lower sensitivity and poorer velocity resolution, the C 18 image is barely resolved 
and the velocity pattern is sufficiently distorted to make analysis difficult. 

In the left panel of Figure 2, contours of 2.7 mm continuum emission at resolution 2'.'2 are 
overlaid on the color image of 13 CO integrated intensity for AB Aur. Line and continuum 
emission peaks are spatially coincident to within 0'/3. The optically thin nature of the 
2.7 mm emission results in a direct relationship of disk mass to flux, Md=l00(F u d 2 / k v B v (T)). 
Following MS97, we calculate the disk mass using this expression and the 11.5 mJy 2.7 mm 
flux, assuming a gas to dust mass ratio of 100 and ^=0.009 cm 2 g _1 at 2.7 mm. Adopting 
a temperature of 40 K, we derive M d =0.009 M Q with a factor of a few uncertainty. S04 
obtained a similar value, M d =0.013 M , with a factor of 7 uncertainty. The center panel 
of Figure 2 displays our 2.7 mm continuum contours overlaid on the near-IR scattered light 
image of F04, based on which they suggest spiral structure. While the scattered light images 
trace small amounts of material, our 2.7 mm image demonstrates that the spiral features are 
massive and represent density contrasts, as the 2.7 mm flux directly traces the mass. The 
three most intense scattered light features are detected at least at the 3a level of significance 
in 2.7 mm continuum, while the strong northern and southeastern scattered light arms are at 
better than the 5a level. These northern and southeastern features are more obvious in the 
right panel of Figure 2. There, the symmetric central core of the continuum image, containing 
more than 75% of the total flux, has been subtracted to enhance the outer regions. Given 
the uncertainties in the symmetric source fit, the resulting northeastern and southeastern 
asymmetry residuals provide between 5% and 11% and 7% to 14% of the total emission, 
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respectively; residual peak locations are uncertain by less than 0M, and emission peaks are 
at least at the 3a level. 

The smooth velocity gradient in the middle panel of Figure 1 suggests that material is 
bound to the star and executes Keplerian orbits, a hypothesis supported by S04. The top 
panel of Figure 3 displays the morphology of the 13 CO emission in velocity intervals of width 
0.34 km s _1 . The resolution of these robustly weighted maps, 2'.'2 xl'.'8, is almost a factor 
of three better than achieved by MS97 and S04, as is the signal-to-noise ratio. Continuum 
resolution and signal-to-noise are similarly improved. Robust weighting allows us to em- 
phasize small scale features in individual velocity bins as opposed to global characteristics. 
Enhanced resolution and sensitivity enable a detailed comparison of the line emission at 
different velocities with that expected from models. The model consists of a flat, Keplerian 
disk with a single component power law emission profile and an outer radius cutoff. The 
emission profile was fit empirically and the emission was binned by line-of-sight velocity and 
convolved with the 2 / . / 2xl'/8 beam. The observed emission was subtracted from the model 
and the goodness of fit was measured with a difference squared merit function. Our best 
fit model and 2a clipped residuals are presented in the middle and bottom panels of Figure 
3. Our interpretation relies entirely on the 13 CO fit. Given the poorer angular and spec- 
tral resolution of the 12 CO and C 18 observations, model fits for those lines were used as 
consistency checks only. The upper section of Table 1 lists the best fit parameters for each 
spectral line. 

The best fit dynamical mass for the 13 CO data is 2.8±0.1 M Q . Disk radius, inclination, 
and PA are 6151^ AU, 21.5^3° and 58.6±0.5° respectively. The emission profile has form 
S = 48.4(r/100)- L28 mJy inside a radius of 615 AU. Results for the 12 CO and C 18 are 
reasonably consistent. We discuss any large discrepancies below. Estimates using other 
methods to obtain i, PA and disk radius are shown in the lower section of Table 1. The 
uncertainties quoted throughout this work are purely statistical estimates, do not include 
possible errors in the distance, and are 3a for our model fits and la for other methods. If 
the velocity pattern is not Keplerian the errors on all disk parameters except radius could 
be substantially larger. Disk radius is more correlated with the assumed emission profile. 

3. Discussion 

Our resolved images of the gas and dust emission from AB Aur enable us to go be- 
yond earlier millimeter observations at lower resolution and determine a central mass of 
2.8±0.1 M . This dynamical mass is only a little higher than the 2.4±0.2 M Q obtained from 
photometry and pre- main sequence track fitting (van den Ancker, de Winter, & Tjin A Djie 
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1998). Our dynamical model results in several local x 2 minima but these are significantly, 
statistically deviant (>20er) from the global minimum value. Our best fit values of i, PA, and 
disk radius (see Table 1) are more accurate than those derived by S04 from lower resolution 
data and more complex models, and well within their uncertainties. 

Large scale envelope emission, present for r>600 AU (S04), may contaminate more 
optically thick CO emission. The best fit 12 CO emission profile described in Table 1 indicates 
structure extending to beyond 1000 AU and flux ratios of CO isotopomers show that 12 CO 
is optically thick. Given this, the high reduced x 2 value (~7) for our 12 CO model fit is 
not surprising and demonstrates that a pure disk model does not fit the 12 CO emission. 
The largest deviation of 13 CO from the model (~4.5<r) occurs in the second velocity bin 
(see Figure 3). The deviation is hardly significant and could be due to a number of factors 
including local density enhancement, a flared disk surface, etc. A detailed discussion of 
these is beyond the scope of this work, but there may be some evidence of disk flaring in 
the top panel of Figure 3 where the southern lobe of emission in bins 4 through 7 splits 
into two pieces. Given that the 13 CO observations fit the pure disk model much better than 
12 CO and there is known contamination in the 12 CO emission, discrepancies between disk 
properties derived from 12 CO and 13 CO can be wholly attributed to this effect, increasing 
our confidence in the 13 CO values. 

The model fits to mass and i based on 13 CO and C 18 emission are consistent within 
the errors and the difference in radii is not unexpected given the relatively low abundance of 
C 18 0. Our PA determinations are consistent with most published values. Inclinations in the 
upper section of Table 1 are consistent with near-IR interferometry results (Millan-Gabet 
et al. 2001; Eisner et al. 2003, 2004) for the inner 1 AU of the disk, and with outer disk 
values from S04. We checked our model values by calculating i using only the separation of 
the extreme velocity emission and, as Table 1 shows, again find good agreement. However, % 
values determined from axis ratios in the visibility plane are systematically higher than the 
velocity based methods. Recent determinations of i=45-65° and PA=30±15° by Liu et al. 
(2004) from mid-IR adaptive optics observations with nulling interferometry, agree less well 
with our model derived values, and their inclination is much closer to the larger, axis-ratio- 
determined values of Table 1. We believe that discrepancies between the values of inclination 
and position angle derived by different methods are probably due to assumptions about AB 
Aur symmetry. 

S04 and F04 demonstrate that the assumption of symmetry is invalid, as do our mil- 
limeter observations. The spiral structure, initially suggested by F04 and demonstrated to 
be a density contrast through this work (see Figure 2), is a major source of asymmetry. If 
spiral structure is present and streaming motions along the arms are significant, the gas dis- 
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tribution could be tightened along the arm, resulting in foreshortening perpendicular to that 
direction. Whether inclination increases or decreases depends on source geometry. In either 
case, the PA will be biased toward the direction of the long axis of the spiral arm. Spiral 
structure could also affect the % and PA values determined from velocity structure, depend- 
ing on the strength of the features and the optical thickness of the line; more optically thin 
emission features would experience greater influence from such local density enhancements. 
Sub-Keplerian rotation will generally result in larger inclinations for velocity determined 
methods with fixed mass, but a simple gas pressure aided sub-Keplerian disk is probably too 
simplistic an assumption if streaming motions are altering the bulk of the gas motion along 
the arms. 

Our observations show two effects that could be due to the influence of spiral structure: 
1) the model fit to the optically thin C 18 line leads to a significantly greater PA than does 
the 13 CO fit. The direction of change of the PA is along the strong southeastern arm as 
would be expected if local streaming velocities were playing a role; 2) in the left panel of 
Figure 2, the 13 CO emission spreads well beyond the boundaries of the continuum emission 
in the northeast and southwest but has similar extent in the southeast. If spiral structure 
is changing gas morphology it would cause larger spatially determined inclinations and PAs 
thereby explaining the higher inclinations estimated from axis ratios. We conclude that 
the range of published values of inclination and PA can be explained by the presence of 
substructure. The values derived from our kinematic models may be influenced by this 
substructure, but we do not expect the effects to be great. 

The wealth of observations of AB Aur at different wavelengths and at different spatial 
and spectral resolutions enables analysis of this complicated circumstellar disk from the 
inner 1 AU (e.g. Eisner et al. 2004) to the outer envelope beyond 1000 AU (e.g. S04). Our 
detailed observations give an accurate dynamical mass and resolve issues about i and PA, 
indicating that apparent inconsistencies are due to asymmetry. Our high resolution images 
in 2.7 mm continuum and 13 CO and C 18 gas emission, show evidence of spiral structure. 
Such structure has important implications for planet formation. If it is due to instability 
and persists long enough or is sufficiently strong, local collapse can occur, quickly forming 
large planets (Boss 2002; Pickett et al. 2003 and references therein). Alternatively, if the 
instability dissipates or the structure is due to existing planets (Bate et al. 2003), the increase 
in local density may increase the cross section for collisions of planetesimals and speed up 
the core accretion process (Rice et al. 2004). 

Research with the OVRO array is funded in part by NSF grant AST-9981546. SC is 
supported by an NSF Graduate Research Fellowship and JE by a Michelson Fellowship. We 
thank M. Fukagawa for providing her scattered light image of AB Aur, J. Carpenter for 
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Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 

Fig. 1. — OVRO millimeter array images of the 12 CO (left panel), 13 CO (center panel) and 
C 18 (right panel) emission from AB Aur. Intensity-averaged velocity maps (first moment), 
shown in color, are overlaid with contours of the integrated intensity (zeroth moment). 
Contours begin at the 3a level and are spaced by 3a for 12 CO and 13 CO. For C 18 maps 
spacings are 2a. The integrated 12 CO, 13 CO, and C 18 maps have la rms of 145, 36, and 
39 mJy km s _1 respectively. Emission was clipped at least at the 4a level in the first moment 
maps. Circularly restored naturally weighted equivalent beam sizes are shown in the bottom 
left corners. 




Right Ascension (J2000) Right Ascension (J2000) Right Ascension (J2000) 



Fig. 2. — Left panel: 2.7 mm continuum emission for AB Aur (contours) overlaid on the 
13 CO integrated intensity image (color). The lower left portion of the map shows the con- 
tinuum beam size. The 13 CO beam surrounds the continuum beam in the left panel. Visible 
13 CO emission begins at 3a level while the continuum contours begin at 3a and increase 
as 2a where la is 0.36 mJy per beam. The center panel shows the same contours of this 
continuum emission on the scattered light image of F04. The right panel shows contours 
of the continuum image with a best-fit gaussian source subtracted. Contours for this panel 
begin at 3a and increase as la. 
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Fig. 3. — Red, yellow, green, blue and dark blue emission represent approximately 2, 3.5, 
5.0, 6.5 and 8.5a above the noise. The left-most images have velocity 4.27 km s _1 . Image 
velocity width is 0.34 km s _1 . Large tick marks define the position a =04 h 55 m 46 s and 
5 =30°33'05". Right ascension ranges from 46.58 s to 45.23 s left to right and the declination 
extends from 32'54.8" to 33'12.3" bottom to top. Upper Panel: 13 CO observations with 2"2 
by 1'.'8 resolution and la rms noise of 35 mJy/beam. Middle Panel: Best fit Keplerian model 
with mass=2.8 M , 2=21.5°, and PA=58.6°. Bottom Panel: Residuals of the data and model 
clipped at 2a. 



Table 1. Best-fit Model Parameters 



Line Mass Radius Inclination Position Angle Method 

Mg AU Degrees Degrees 



13 co 
c 18 o 
12 co 


2.8±0.1 

„ 7y +0.09 
z - ' ' -0.15 

3 25 +o.i3 


615+* 
49711° 
1060±10 


21 5+°' 4 
21 4+°' 7 

on tr+0.4 

32.5_ 6 


58.6±0.5 
75±2 
61±3 


Model 
Model 
Model 


13 co 


2.79 


N/A 


24±4 


69±8 


Extreme Velocity 


c 18 o 


2.79 


N/A 


17±6 


73±21 


Extreme Velocity 


12 co 


2.79 


N/A 


26±9 


66±25 


Extreme Velocity 


13 co 


N/A 


390±15 


33±5 


85±10 


UV Plane 


12 co 


N/A 


500±13 


41±3 


26±4 


UV Plane 


Dust 


N/A 


170±17 


44±12 


79±19 


UV Plane 



Note. — Dust is the 2.7 mm continuum emission. Error bars are 1<t for all but the 
line model where errors are 3<r. Radii are HWHM values except for the line models 
where they represent outer edges of the disk. Emission profile for the line model was of 
the form S=S ( W £ AU )~ k mjy for r<R ou t. Values of (S ,k) are (48.4±J;f ,1.28±0.04), 
(13.4l^; 2 ,1.30lg;°®), and (150±2,1.1±0.1) for 13 CO, C ls O, and 12 CO respectively. 



